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The issues involved in a determination of the baryon resonance spectrum in lattice QCD are discussed. The 
variational method is introduced and the need to construct a sufRcient basis of interpolating operators is em- 
phasised. The construction of baryon operators using group-theory techniques is outlined. We find that the use 
both of quark-field smearing and link-field smearing in the operators is essential firstly to reduce the coupling of 
operators to high-frequency modes and secondly to reduce the gauge-field fiuctuations in correlators. We conclude 
with a status report of our current investigation of baryon spectroscopy. 



1. INTRODUCTION 

Spectroscopy is a powerful tool for uncover- 
ing the important degrees of freedom of a phys- 
ical system and the interaction forces between 
them. The spectrum of QCD is very rich: con- 
ventional baryons (nuclcons, A, A, S, f2, etc.) 
and mesons (tt, K, p, etc.) have been known 
for nearly half a century, but other, higher-lying 
exotic states, such as glueballs, hybrid mesons 
and hybrid baryons bound by an excited gluon 
field, and 'multi-quark' states, consisting predom- 
inantly of four or five quarks in the case of mesons 
and baryons respectively, have proved more elu- 
sive, partly because our theoretical understanding 
of such states is insufñcient, making their identi- 
fication difficult. 

Interest in excited baryon resonances in par- 
ticular has been sparked by experiments dedi- 
cated to mapping out the N* spectrum in Hall 
B at the Thomas Jefferson National Accelera- 
tor FaciUty (JLab). Much of our current under- 
standing of conventional and excited hadron res- 



onances comes from QCD-inspired phenomeno- 
logical models. For conventional baryons, the 
extensive calculations by Isgur, Karl, and Cap- 
stick within a non-relativistic quark model |fTl2l3| 
remain inñuential. However, there are a growing 
number of resonances which cannot be easily ac- 
commodated within quark models. States bound 
by an excited gluon field, such as hybrid mesons 
and baryons, are stiU poorly understood. The na- 
tures of the Roper resonance and the anomalously 
light A(1405)~ remain controversial. Experiment 
shows that the first excited positive-parity spin- 
1/2 baryon lies below the lowest-lying negative- 
parity spin-1/2 resonance, a fact which is difficult 
to reconcile in quark models. The question of the 
so-called "missing" baryon resonances is still un- 
resolved: the quark model predicts many more 
states0E] than are currently known. Compared 
to the large number of positive-parity states, 
there are only a few low-lying negative-parity res- 
onances. A quark-diquark picture of baryons pre- 
dicts a sparser spectrum^]. 

Given the current intense experimental efforts 
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in spectroscopy for baryons, the need to predict 
and understand the baryon spectrum from first 
principles is clear; lattice QCD calculations pro- 
vide the means of undertaking such ab initio stud- 
ies. The aim is not merely to obtain a set of 
masses for the states, but also to gain insight into 
the quark and gluon structure of the states and to 
understand the relevant degrees of freedom; this 
latter aspect will be an important emphasis of this 
talk. The freedom to vary quark masses, numbers 
of quark flavors and even the gauge group enables 
us not only to relate lattice computations directly 
to experiment, but also to QCD-inspired pictures 
of baryon structure. 

The layout of the remainder of this talk is as 
foUows. In the next section we introduce the vari- 
ational method as a means of extracting infor- 
mation about the QCD spectrum, and demon- 
strate the need to construct a sufficient basis of 
interpolating operators. Section |31 describes the 
construction of lattice baryon operators using a 
group-theory methodl^l; an alternative approach 
we havc dcvcloped is described in ref. W. 

2. HIGHER EXCITED RESONANCES 
AND THE VARIATIONAL METHOD 

A comprehensive picture of resonances requires 
that we go beyond a knowledge of the ground 
state mass in each channel, and obtain the masses 
of the lowest few states of a given quantum num- 
ber. This we can accomphsh through the use of 
the variational method'7J8j. Rather than mea- 
suring a single correlator C(í), we determine a 
matrix of correlators 

where {Oí; i = 1, . . . , N} are a basis of interpolat- 
ing operators with given quantum numbers. We 
then solve the generahzed eigenvalue equation 

C{t)u ^ X(t,to)C{to)u 

to obtain a set of real (ordered) eigenvalues 

^n{t,to), where Aq > Ai >, At large Eu- 

chdean times, these eigenvalues then dehneate be- 
tween the different masses 

A„(¿,ío) ^ g-AÍ„(í-í„) ^Q(g-M„+i(t-t„))^ 



The eigenvectors u are orthogonal with metric 
C(ío), and a knowledge of the eigenvectors can 
yield information about the partonic structure of 
the states. A recent application of this method to 
the baryon spectrum using operators constructed 
using smeared-quark sources of varying widths is 
discussed in ref. 

3. BARYON OPERATORS AND THE 
LATTICE 

Crucial to the apphcation of variational tech- 
niques is the construction of a basis of opcrators 
that have a good overlap with the lowest-lying 
states of interest. These operators should have 
the property that they respect the symmetries of 
the lattice, rather than being mere discretisations 
of continuum interpolating operators. The LHP 
CoUaboration has developed techniques to enable 
the construction of baryon interpolating opera- 
tors that can easily be extended to include multi- 
quark operators, and those with excited glue 5.6j . 

3.1. Interpolating operators and lattice 
symmetries 

States at rest are classified according to their 
transformation properties under rotations; in a 
lattice calculation, such rotations are restricted 
to those of the cubic group of the lattice, O. This 
group has the foUowing properties: 

• O has 24 elements 

• There are five conjugacy classes, and 
hence five single-valued representations, 
Ai, A2, E,Ti, and T2, of dimensions 1,1,2,3 
and 3, respectively. 

The irrcducible representations are defined such 
that, under the elements R of O, the operators 
lying in thc A irreducible representation (irrep.) 
transform as 

L/(i?)0fV(i?)t = ^0Í^)l?i^(ñ) 

A' 

where A is the row of the irrep., U is the uni- 
tary matrix effecting the rotation, and D is the 
corresponding representation matrix. 

The addition of the spatial inversion opera- 
tor, corresponding to parity, yields the Octahedral 
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Table 1 

Continuum limit spin identification: the number nj^ of times that the A irrep. of the octahedral point 
group Oh occurs in the (reducible) subduction of the J irrep. of SU{2). The numbers for GiujG^u, 
are the same as for Gig, G^g, Hg, respectively. 
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Group Oh, and the irreps. acquire an additional 
subscript g or w, denoting positive and negative 
parity respectively. In the case of baryons, we 
must consider the double-valued, or spinorial, ir- 
reducible representations. There are three such 
irreducible representations, Gi, G^ and H, of di- 
mensions 2, 2 and 4 respectively, again with g and 
u labels to denote parity. 

The irreducible representations J of the con- 
tinuum rotation group SU(2) are reducible under 
the cubic group O; the number of times that 
the A irrep. of Oh occurs in the reducible sub- 
duction of thc J irrep. of the continuum SU(2) 
is shown in Tablenj States with J > 5/2 He in 
irreducible representations containing states with 
lower spins, and furthermore, for a given J, differ- 
ent continuum hehcities can correspond to mem- 
bers of different irreducible representations of O/,,. 
The masses of the components in these distinct 
irreps. wiU agree only in the continuum hmit. 

An imphcit assumption in previous lattice 
studies is the increase in ground-state masses with 
increasing spin. This assumption is not necessar- 
ily reahzed in nature; in the nucleon sector, the 
lowest-lying = f ^ state is comparable in mass 

to the lowest-lying J^ = |^ state^üj. In a lattice 
computation, the spin of an energy level in a given 
channel can only be identified by an examination 
of the degeneracies between energies in different 
irreps. in the approach to the continuum hmit. 
Thus the ability to extract several energy levels 
in each channel, and therefore the work outhned 
in this talk, is especiaUy cruciaL 



3.2. Operator recipe 

The starting point for the construction of our 
three-quark operators is a basis of gauge-invariant 
terms of the form 

$;^,^4-.=^-..c(A^^^^)l(¿f^)f,(¿i''^^)f,, (1) 

where A, B, C indicate quark ñavor, a, b, c are 
color indices, a,P,"f are Dirac spin indices, ip in- 
dicates a smeared quark field, and ¿j*'-' denotes 
the p-hnk covariant displacement operator in the 
j-th direction; the quark fields are smeared using 
a three-dimensional gauge-covariant Laplacian. 

The displacements in eqn. are chosen so 
as to span a range of possible spatial configura- 
tions of the quarks, and are listed in Table|21 In 
particular, the singly-displaced operators mimic 
a quark-diquark combination, whilst the doubly- 
and triply-displaced operators are chosen since 
they may favour A- and K-ñux configurations, 
respectively. These operators are designed to al- 
low a large number of baryon operators to be con- 
structed using only a small number of quark prop- 
agators, by exploiting the cubic symmetries of the 
correlators to reduce the number of quark sources 
that need to be employed. 

These gauge-invariant operators are now com- 
bined into elemental operators having the appro- 
priate flavor structure; we assume exact isospin 
symmetry, and classify operators according to 
their isospin, and strangeness (or charm etc. 
where appropriate) . We thus obtain a set of 
gauge- and translational-invariant elemental op- 
erators Bf{t) — J2x i^'^) having the appro- 
priate flavor structure. The flnal step in our pro- 
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Table 2 

The six types of three-quark ^aPr, ijk opera- 
tors used, where A,B,C indicate the quark fla- 
vors, 1 < a,/3,j < 4 are Dirac spin indices, 
and —3 < i,j,k < 3 are displacement indices. 
Smeared quark fields are shown by soUd circles, 
hne segments indicate covariant displacements, 
and each hollow circle indicates the location of a 
color Eabc couphng. A displacement index having 
a zero value indicates no displacement. 

Operator type Displacement indices 



i — j — k — 

í = j = 0, fc ^ 

¿ = 0, j = -fc, /c ^ 

í = 0, IjI y¿ \k\, Jk ^ 

í = -J, \3\^\kV jk^O 
m^\j\^\k\, ijk^O 



single-site 

singly-displaced 
doubly-displaced-I 



-displae 



doubly-displaced-L 



triply-displaced-T 
triply-displaced-0 



cedure is to apply the group theoretical projection 
to yield a set of operators B^'^^ that transform 
according to thc row A of the A irreducible rep- 
resentation: 



B 



AXF 



^•^h Reol 



r[''^{R)UnB[ml (2) 



where Of^ is the double group of Oh, R denotes an 
element of O^ , gQO is the number of elements in 

O^ , and dj^ is the dimension of the A irreducible 
representation. 



4. IMPLEMENTATION 



sive study of baryon spectroscopy in lattice QCD. 
We begin by describing the tuning of the smear- 
ing Darameters¡ll|. 

4.1. Smearing procedure 

It has long been known that the damping of the 
couphngs to short-wavelength, high-frequency 
modes is a crucial requirement for the extraction 
of the masses of the lowest-lying states in a lattice 
QCD calculation j 1 21 1 3| . In our study, we employ 
a gauge-covariant Gaussian smearing: 



where 

(Uk{x)^{x + k)-^{x)y (3) 

fc=±l,±2,±3 

There are two tunable parameters, the smearing 
radius cts a-nd the number of iterations n^. 

The operators introduccd in cqn. ^ in gen- 
eral involve gauge-covariant displacements. In or- 
der to reduce the statistical fluctuations in the 
gauge fields, and possibly to further damp out 
the high-frequency couphngs, once can smear 
the hnk variables, both in the Laplacian and in 
the displacement operators. A widely adopted 
procedure for so doing is that proposed by the 
APE coUaboration 14 , which involves a projec- 
tion back onto SU{3). We instead adopt the an- 
alytic "stout" smearing prescription proposed in 
ref. ¡15¡, defincd by the iterative procedure 

Ut+'\x) = exp(ípe(")(a;)){/i")(x), (4) 

Qk{x) = '-{ni{x)-nk{x)) - 

^TT(nl{x)~nk{x)) (5) 

nk{x) = Ck{x)Ul{x) (6) 
Ck{x) = ^ (c/,(:E)C/fe(a; + í)C//(x-|-fc) 

+U¡{x - Í)Uk{x - Í)U,{x - i + k)) . (7) 



The remainder of this talk wiU provide a status 
report of our program to implement a comprehen- 



Here there are again two tunable parameters, the 
number of iterations iip, and the staple weight p. 
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Figure 1. EfFective masses M{t) for unsmeared (circles) and smeared (triangles) operators 
0s5, OsD, Otdt- Top row: only quark-field smearing n„ = 32, Us — 4.0 is used. Middle row: 
only Hnk-variable smearing Up = 16, Upp = 2.5 is applied. Bottom row: both quark and link smearing 
rirj = 32, as = 4.0, Up = 16, Upp = 2.5 are used. Results are obtained on 50 quenched configurations on 
a 12^^ X 48 anisotropic lattice using the Wilson action with a¿. ^ 0.1 fm, as/at = 3.0^2], as described in 
the text. 
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For our tests of the efficacy of quark-field and 
gauge-fink smearing, we employed 50 configura- 
tions of a quenched, anisotropic 12^^ x 48 lattice, 
using the standard Wilson gauge and fermion ac- 
tions. The spatial lattice spacing, determined 
from the string tension, is ag ~ 0.1 fm, with a 
renormahzed anisotropy as/at = 3; the quark 
mass was chosen so that rn^ ~ 700 MeV. Correla- 
tors were computed corresponding to the single- 
site, singly-displaced and triply-displaced-T op- 
erators of Table [21 for the single-site operator, 
a projection onto Gig was performed, whilst for 
the remaining operators a single Dirac componcnt 
was chosen with no attempt to project onto an ir- 
reducible representation. 

The quark-field smearing parameters ria- and a 
were determined by requiring that the effective 
mass for these operators reach a plateau as close 
to the source as possible. The gauge-hnk smear- 
ing parameters were tuned so as to minimize the 
noise in the effective masses; for these lattices we 
found optimal smearing parameters n^p = 2.5 
with Up = 16. The use of gauge-link smearing 
has only a smaU effect on the mean values of the 
baryon effective masses. However, we found a 
dramatic reduction in the statistical variance in 
the singly-displaced and triply-displaced opera- 
tors, as demonstrated in Figuren 

4.2. Group theory projections 

Having demonstrated the efficacy both of hnk- 
and quark-smearing in isolating the ground-state 
energies in the correlators, we now proceed to 
apply the projection formula of eqn. Ef- 
fective masses obtained from three selected nu- 
cleon operators in the Gig, Giu and Hg channels, 
computed using the singie-site, doubly-displaced- 
I and triply-displaced-T elemental operators, are 
shown in Figurc[21 

An cxploratory study of the nucleon spectrum 
using the Clebsch-Gordon methodp] was per- 
formed in ref. J7]. Here a larger 16'^ x 64 lattice 
was employed using around 300 configurations 
with a quark mass — 500 MeV. Such an en- 
semble was sufficient to enable the appfication of 
the variational method to extract the lowest-lying 
eigenvalues of the Gig correlator matrix, as shown 
as Figure|31 Furthcrmore, a plateau was found in 




Figure 3. The lowest-lying effective masses in 
the Gig channel obtained using the variational 
method[T7]. 



the G^g irrep. one accessible only through the 
use of displaced, rather than single-site, opera- 
tors. The effective masses of the ground states in 
each of the positive-parity Gig,Hg and G^g irre- 
ducible representations are shown in Figure01 the 
apparent coincidence between the Hg and G^g ef- 
fective masses suggests that the identification of 
the lowest-lying Hg nucleon state with spin-3/2 
is somewhat premature, and enforces the need for 
the program outhned in this talk. 

5. CONCLUSIONS 

We have outhned a program to study the res- 
onance spectrum in lattice QCD. The use of the 
variational method and the need to isolate several 
energy levels in each channel require a sufficiently 
broad basis of operators. Having developed suit- 
able group-theory mcthods to project operators 
onto the irreducible representations of the cubic 
group, and having examined the efficacy of both 
quark- and gauge-hnk-smearing, we are now iden- 
tifying a more hmited set of operators that we will 
employ in a large-scale study of the hadron spec- 
trum. Our methods are apphcable not only to 
baryons, but also to mcsons, to states with ex- 
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Figure 2. EfFective masses for three selected nucleon operators: a single-site operator in tlic 1=1/2 
Gig cliannel (left), a doubly-displaced-I operator in the Gi„ channel (center), and a triply-displaced-T 
operator in the Hg channel. The smearing parameters used were no- = 32, = 4.0, Up = 16, Upp = 2.5. 
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Figure 4. The ground-state effective masses in 
the Gig, G^g and Hg irreducible representations. 



cited glue, and to multi-quarlc and multi-hadron 
states. Only by performing such a program can 
we hope to identify the states of QCD, and in par- 
ticular their spins and parities, in the continuum 
limit. 
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